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ABSTRACT 

The supersymmetric QCD corrections to the hadronic width of a heavy charged 
Higgs boson, basically dominated by the top-quark decay mode —^tb, are 
evaluated at 0{as) within the MSSM and compared with the standard QCD 
corrections. The study of such quantum effects, which turn out to be rather 
large, is essential to understand the hypothetical supersymmetric nature of 
a heavy charged Higgs boson potentially produced in the near future at the 
Tevatron and/or at the LHC. 



After the discovery of the top quark at the Tevatron |I| , the "next-to-leading" task in 
order of importance in Particle Physics is the search for Higgs bosons. Indeed, whereas 
the finding of the top quark completes the spectrum of matter fields in the Standard 
Model (SM) and it was a necessary condition for the SM to be correct, it is not at all a 
sufficient condition. The decisive proof lies in the nature of the spontaneous symmetry 
breaking mechanism (SSB). Thus the very crucial question is still left pending: is the 
SSB caused by truly fundamental scalars or is it triggered by a dynamical mechanism 
involving a new species of fermions and/or a new strong interaction force?. This question 
is as old as the proposal for the first serious extensions and alternatives to the SM, namely, 
Supersymmetry (SUSY) and Technicolour. Here we will be concerned with the Minimal 
Supersymmetric Standard Model (MSSM) 0. 

In the near and middle future, with the upgrade of the Tevatron, the advent of the 
LHC, and the possible construction of an e~ supercollider, new results on top quark 
physics, and perhaps also on Higgs physics, will be obtained in interplay with Supersym- 
metry that may be extremely helpful to complement the precious information already 
collected at LEP from Z-boson physics. Here we wish to dwell on the phenomenology 
of supersymmetric Higgs boson decays into hadrons with an eye on these future develop- 
ments. Two basic parameters are needed to deal with the (tree-level) MSSM Higgs sector, 
namely (M^o, tan/5) |^. In the present study (see also the companion paper @]), it will be 
useful to divide the supersymmetric Higgs parameter space into two mass segments: light 
CP-odd higgses (M^o < Mz) and heavy CP-odd higgses (M^o > Mz). However, since 
these segments are unequally sensitive to large and small values of tan (3, it will prove con- 
venient to define the following two relevant regions of the (M^o, tan/?)-parameter space. 
One of them, call it Region I, is characterized by Higgs masses in the light segment to- 
gether with large values of tan/3 of order ~ rrit/mh. The other region (Region II) consists 
of Higgs masses in the heavy segment in combination with more moderate values of tan (3, 
say m the interval 2 < tan/5 < 30[|. 

In Refs. 1^, whose notation and definitions we shall adopt hereafter, we studied 
the impact of Region I on the top-quark decay process. Here we will focus on a related 
process, namely, the top quark decay mode of a charged supersymmetric Higgs, ^ tb, 
which is allowed in Region II, provided that M^o > 160 GeV. We shall again be partially 
concerned with Regions I and II in the companion paper 0] , where we round off our study 
by addressing the quark-antiquark decay modes of the three neutral Higgs states of the 
MSSM. 

^ Incidentally, we note that Regions I and II correspond (approximately) to domains of the 
{Mao , tan /3)-plane formerly used [|| to alleviate some of the -still not completely banished- Z-boson 
anomalies, where in addition it was required the concurrence of a chargino and a stop both as light as 
permitted by phcnomenological bounds, namely, of about 60 — 70 GeV. 
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While a simple tree-level study of ^ tb is blind to the nature of the Higgs sector 
to which belongs, it is clear that a careful study of the quantum effects on iJ+ — > 1 6 
could be the clue to unravel the potential supersymmetric nature of the charged Higgs; in 
particular, it should be useful to distinguish it from a charged Higgs belonging to a general 
two-Higgs doublet model. The charged Higgs boson can decay hadronically into several 
quark final states, and if it is sufficiently heavy it can also decay into top and bottom 
quarks. It is natural to tackle the radiative effects on ^ tb hj first considering the 
QCD corrections. The conventional QCD corrections were first considered in Ref. [§] and 
they are known to be large and negative in the limit of Higgs masses much bigger than the 
quark masses. In spite of the huge size of the standard QCD corrections for light quark 
final states, it is clear that these channels are severely suppressed by the small Yukawa 
couplings grriq/Mw <^ l(forg = u,d,c,s) and/or by off-diagonal CKM matrix elements, 
so that their branching fractions are very tiny. Therefore, as soon as the tb threshold 
is open {Mh+ > rrit + mi, ~ 180 GeV) one is left with tb as the only relevant 

hadronic decay of a heavy charged Higgs boson. The conventional QCD corrections 
to that decay ^ cannot distinguish the nature of the underlying Higgs model, but 
their knowledge is indispensable to probe the existence of additional sources of strong 
quantum effects beyond the SM. Here we will concentrate on the SUSY-QCD quantum 
effects mediated by squarks and gluinos and shall compare them with the standard QCD 
corrections. A few remarks on the remaining MSSM corrections will be made later on. 

The relevant SUSY-QCD one-loop decay diagrams constructed from gluinos and squarks 
(stop and sbottom species) are seen in Fig. la. In Fig. lb we sketch one possible mechanism 
that could be used at the LHC to annually produce ~ 10"^ — 10^ charged Higgs particles 
of a mass comprised between a few hundred GeV up to about 1 TeV for a luminosity of 
C ~ 10^^ cm~^ |T^. There are both the tb- and 6t-fusion mechanisms, which are to 



be treated together with the processes gb —>■ i and gb —>■ H~ t as carefully explained 
in|]lT|. Trigging on a top quark in association is very useful to avoid the signal being 
swamped by huge backgrounds. These mechanisms are primarily initiated by two highly 
energetic gluons, and use as a production vertex the very same decay vertex that we wish 
to study. We point out, in passing, that these diagrams also contribute to the cross-section 
for single top-quark production, whose measurement is one of the main goals at the next 
Tevatron run (Run II). It is not our aim to give the complete list of diagrams, but it is 
clear that these production mechanisms could be rather efficient in the colliders, for the 
1 6- vertex can be strongly enhanced and result in a very distinctive phenomenology 
as compared to the experimental expectations for the SM Higgs production, typically 
through (one-loop) g g-insion [|TU|. Most important, in the MSSM the if"*" 1 6- vertex can 
receive significant corrections (in some cases of order 50%) which could play a decisive 
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role to disentangle whether a charged Higgs hypothetically produced in a hadron collider 
is supersymmetric or not. 

To evaluate the relevant corrections to F = r{H^ —>■ tb) in the MSSM, we shall 



adopt the on-shell renormalization scheme 12] where the fine structure constant, a, and 



the masses of the gauge bosons, fermions and scalars are the renormalized parame- 
ters: (a, Mw, Mz, MH,mf, Msusy, •••)• The interaction Lagrangian describing the tb- 
vertex in the MSSM reads as follows: 

'^Htb = l^!' H"^ i [rrit cot (3 Pl + ruh tan /3 Pr] 6 + h.c. , (1) 
V 2MiY 

where Pl,r = 1/2(1 =F 75) are the chiral projector operators, tan (3 is the ratio between 
the vacuum expectation values of the two Higgs doublets of the MSSM and Vtb is the 
corresponding CKM matrix element-henceforth we set Vtb = 1- 

The on-shell renormalized vertex tb in Fig. la is derived following a straightforward 



generalization of standard procedures in the SM |jT2|] and can be parametrized in terms of 
two form factors Hi, Hr and the corresponding mass and wave-function renormalization 
counterterms 5raf and 5Z{^ ^ associated to the external quark lines in Fig. la: 



with 



iG= [m* cot /3(l + Gi.)Pz. + mfe tan/5 (l + G'H)Pij] , (2) 

V 2 My/ 



The counterterms are computed in the on-shell scheme from the SUSY-QCD interaction 
Lagrangian involving squarks and gluinos 0. The explicit results, which we refrain from 
repeating here, are common to those in Ref. |]^. As for the vertex form factors, they 
depend not only on pure strong supersymmetric interactions but also on the Higgs-stop- 
sbottom semiweak interaction Lagrangian [Q. One finds (summation is understood over 
indices a,b = 1, 2): 

Tflf CO u l) 

rrib tan p 

(4) 

where C... = C ,__{P-,P' )^gi'^ti,-i^bc) standard three-point functions also carrying 
indices a, 6 summed over; Cp = {N^ — l)/2Nc = 4/3 is a colour factor obtained after 
summation over colour indices, and 

Gab = Rla Rib 9LL + R2a -^26 9RR + -^20 ^Ib 9LR + Rla -^2?) 9RL , (5) 
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with 

qll = sin 2/? — (rrij cot /? + tan/3) , 

gRR = -mjmb(tan/3 + cot/5) , 

9lr = -mb{iJ + Ab tan P) , 

Qrl = -mt{f^ + At cot P) . (6) 
The 2x2 rotation matrices 

i^(^)=f'"1^ 'Z'"] (g = M), (7) 

y sm 9g cos Oq J ^ ^ ^ 

plaguing the previous formulae relate the weak-eigenstate squarks q'^ = {ql, Qr} to the 
mass-eigenstates qa = {qi,q2} as follows: g'^ = J^b^libQb- Therefore the R^"^^ diagonalize 
the corresponding stop and sbottom mass matrices, whose standard form is well-known 
but we will quote here explicitly for convenience: 

A.2_ / ^^+^? + cos 2/3(i- I s2^,)M| mtMij, \ 

^'~[ mtMl^ M4 + m? + |cos2/3.2^M| J ' 

/ + cos2/?(-i + I s2^)M| rribMlj, \ 

' [ nibMl^ Ml + ml- I COS2P Ml J ^ 

with 

MIj^ = At- fi cot f3 , MIj^ = Ab- fi tan p . (10) 

From the renormalized amplitude @ the width of ^ tb, including the one- loop 
SUSY-QCD corrections, is the following: 

r = To {l + ^ [2Re{GL)] + ^ [2Re{Gn)] + ^ [2Re{GL + Gn)]] , (11) 



where the lowest-order result is 

To-f^l^AV^d,^,^), (12) 



with A(l, x^, y"^) = [1 — {x + y)'^] [1 — {x — y)^], and 

D = {Mh+ -ml -ml) [ml cot^ 13 + mltan^ 13) - Amlml, 

Ul = {Ml+ - ml - ml) ml cot^ p , 

Ur = {Ml+-ml-ml)mltan^ p, 

Ulr = -2mlml. (13) 

The numerical analysis of the strong supersymmetric corrections to r{H~^ ^ tb) is pre- 
sented in Figs. 2-5. In all figures where the Higgs mass is fixed, we take Mh+ = 250 GeV, 
and we define = as{Mjj+) by means of the (one-loop) expression 

"'^"'"^^ = (33-2n,)log(M„./A,.,) ' 
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normalized as as{Mz) — 0.12, where nj is the number of quark flavors with threshold 
below the Higgs boson mass Mh+- 

We treat the sbottom mass matrix (^) assuming a non-zero mixing parameter M^^. 
Once fi and tan (3 are given, we use Ai, as one of the inputs. However, for the sake of 
simplicity we assume that 6b = 7r/4, so that the two diagonal entries are equal to a 
common value. We denote by m^^ the lightest sbottom mass-eigenvalue and take it as 
the remaining input; typically, it is bound to satisfy m^^ > 100 — 150 GeV^ from collider 
data, but the limits are not so stringent as those from LEP, namely m^^ > 65GeV. As 
for the stop mass matrix, we choose our inputs as follows. In this case the lightest stop 
mass m^^ is strictly limited only by the LEP 1.5 bound m^^ > 65GeV, and we take it as 
one of the inputs. Thus, upon fixing /x, tan (3 and taking into account that SU (2) ^.-gauge 
invariance requires Mf^ = M^^, it follows that the stop mass matrix depends on one 
remaining parameter which can be taken as At. 

To start with the numerical analysis, we study the dependence of the SUSY-QCD 
effects on the crucial parameter tan f3. In Fig. 2a we plot the SUSY-QCD corrected width, 
eq.(|llD, versus tan/5, for fixed values of the other parameters. For completeness, we 
have included in this figure the partial widths of the alternative decays — >■ Ur and 
— > iy+ /i^, which are obviously free of 0{as) QCD corrections. (To avoid cluttering, 
we have not included — > cs; it is overwhelmed by the r-lepton mode as soon as 
tan (3 > 2.) It is patent from Fig. 2a that, for charged Higgs masses above the t b threshold, 
the decay ^ tb is dominant. Only for very large tan/3 (> 30) and for sufficiently 
big and positive /i (/i > 100 GeV^) -hence outside Region II - the negative corrections to 
—* tb are huge enough to drive its partial width down to the level of r+ z/^. 

Therefore, within the limits of Region II, the top quark decay of the charged Higgs is, by 
far, the most relevant decay mode to look at. 

In considering the various parameter dependences, we present the results of our anal- 
ysis in terms of the quantity 

^, = ^, (15) 
J- 

which gives the relative correction with respect to the tree-level width. At large tan/3, 
the role played by the bottom quark mass becomes very important. Indeed, in Fig. 2b we 
confirm that the external self-energies (basically the one from the 6-line) give the bulk 
of the corrections displayed in Fig. 2a, whereas the (finite) vertex effect is comparatively 
much smaller and its yield becomes rapidly saturated. The existence of potentially large 
SUSY-QCD corrections in the decay tb could be foreseen from the work of Ref. , 

where it is shown that the SUSY-QCD bottom mass corrections are proportional to fi tan /3 
(for fixed, nonvanishing, fi and large tan /?). In our case, these corrections are fed into the 
counterterm Snib/mb on eq.(^ and, when viewed in terms of diagrams of the electroweak- 
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eigenstate basis, they appear as finite contributions proportional to M]^^ generated from 
squark-gluino loop^. 

From Fig. 3a we read off tfie incidence of tlie parameter fi on 6g for various tan (3. We 
see tliat tlie SUSY-QCD correction is extremely sensitive to fi both on its value and on 
its sign. For this reason we have explored a moderate range of fi values. It turns out 
that the sign of the SUSY-QCD correction is basically opposite to the sign of fi, and the 
respective corrections for +/i and for — /i take on approximately the same absolute value. 

Worth noticing is also the dependence on the gluino mass (Fig.3b). The various 
steep falls in that figure are associated to the presence of threshold effects occurring at 
points satisfying rrig + m^^ ~ m^. An analogous situation was observed in Ref. |0 for the 
SUSY corrections to the standard top-quark decay. Away from the threshold points, the 
behaviour of 6g is smooth and perfectly consistent with perturbation theory. On the other 
hand, the sensitivity on m^^ is not dramatic, as can also be appraised in Fig. 3b. This is 
to be expected from the fact mentioned above that it is the bottom (not the top) self- 
energy that dominates the corrections. Virtually for any m^^ , there emerges an important 
correction which raises a long way with the gluino mass before bending -very gently - 
into the decoupling regime, as we have checked. The fact that the decoupling rate of the 
gluinos appears to be so slow has an obvious phenomenological interest. 

In Fig. 4a we analyze in detail the SUSY-QCD correction as a function of the sbottom 
masses. We find convenient to plot 6g versus m^^ , for fixed values of the other parameters 
and for various tan j3. We see that even for m^^ exceeding 250 GeV, and tan (3 > 10, the 
correction remains sizeable {6g > 10%). As another feature, in Fig. 4b we realize that 
Sg is not symmetric with respect to the sign of Af,. Once the sign /i < is chosen, the 
correction is larger for negative values of Af, than for positive values. We have erred on 
the conservative side by choosing Ah = +300 GeV wherever this parameter is fixed. 

In Fig.5 we compare the one-loop SUSY-QCD corrections with the standard 0{as) 
QCD corrections as a function of Mfj+ and different values of tan/5. For the latter 
corrections we use the full analytical formulae of Ref. 00. In the limit Mh+ ^ fnt, the 
standard QCD correction boils down to the simple expression 

r ^ Tqcd - Fq ^ /CV«s\ rnl cot' (| - 6 log ^) + ml tan^ /3 (| - 6 log ^) 
^ Fo V 2 vr y cot^ (3 + ml tan^ f3 

(16) 

This formula is very convenient to understand the asymptotic behaviour. However, as 
we have checked, it is inaccurate for the present range of values of mt unless Mfj+ is 

^In the absence of sbottom mixing, i.e. M^j^ — 0, the large contribution Smi,/mi, ^ — /itan/3 is no 
longer possible, but then the vertex correction does precisely inherits this behaviour and compensates 
for it. However, the condition M^j^ = 0, combined with a large value of tan/3, leads to an scenario 
characterized by a value of Ai, which overshoots the natural range expected for this parameter fl^ . 

^We have corrected several misprints on eq.(5.2) of Rcf. 
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extremely large (beyond 1 TeV) . 

Remarkably, we recognize from Fig. 5 that the supersymmetric QCD effects {6g) can be 
comparable or even larger than the conventional QCD corrections {Sg). For a given tan (3, 
the relative size of the SUSY-QCD effects versus the standard QCD effects depends on 
the value of Mh+- Notwithstanding, it is clear that 6g remains fairly insensitive to Mh+. 

Some discussion on previous work is in order. A first study of the SUSY-QCD correc- 
tions to the hadronic width of a charged Higgs boson is performed in Refs. [|T5[] . However, 
these references either use a too simplified set of assumptions on the spectrum of gluinos 
and squarks (e.g. all squarks are assumed unmixed and degenerate) and/or mi, = is 
assumed, so that the local behaviours of 5g and the crucial effect from finite bottom mass 
corrections at large tan (3 are fully unnoticed. 

Although we have concentrated on the computation of the QCD corrections, it is le- 
gitimate to worry about the larger and far more complex body of electroweak quantum 
effects, especially those coming from possible enhanced Yukawa couplings. A full analysis 
is under way and will be presented elsewhere [|I6|, but a few comments may be necessary 
here. In fact, Yukawa couplings can also give contributions to 6mi,/mi, proportional to 
tan/3 |T3|, but their effect is in general smaller than in the SUSY-QCD case. For ex- 
ample, if there is no large hierarchy between the sparticle masses, the ratio between the 
SUSY-QCD and the Yukawa coupling contributions to Smfy/mfy can be estimated (at high 



tan/?) as Arrig/At times a slowly varying function of the masses of order 1 |13|, where the 
(approximate) proportionality to the gluino mass reflects the aforementioned very slowly 
decoupling rate of the latter. In view of the present bounds on the gluino mass, and 
since At (as well as Af,) cannot increase arbitrarily, we expect that the SUSY-QCD effects 
can be dominant, and even overwhelming for sufficiently heavy gluinos. Notwithstand- 
ing, the sole estimation in terms of Smh/mb may be insufficient, for there are also plenty 
of additional electroweak vertex contributions both from the Higgs sector and from the 
stop-sbottom/chargino-neutralino sector where those Yukawa couplings are also involved. 
Hence, in contrast to the SUSY-QCD case, it is not obvious a priori what is the net 
outcome of the leading electroweak contributions. A partial evaluation of the Yukawa 
coupling effects from supersymmetric Higgs bosons has been made in Ref. and entail 
only a few percent change in the partial width. Even though a complete calculation of 
the Higgs effects - not to mention those associated to the full plethora of supersymme- 
tric particles - is not yet available in the literature, a preliminary analysis made by the 
authors [1^ leads to the conclusion that the over-all leading electroweak supersymmetric 
effects, though not negligible, remain comparatively small and do not drastically alter the 
SUSY-QCD picture presented here. 

In conclusion, we have presented a fairly complete treatment of the SUSY-QCD cor- 



8 



rections to the partial width of the top quark decay of a charged Higgs boson and have put 
forward plenty of evidence that they can be rather large (typically between 10% — 50%), 
slowly decoupling and of both signs. Consequently, they can either reinforce the conven- 
tional QCD corrections or counterbalance them, and even reverse their sign; the QCD 
corrections would then be found much "larger", "missing" or with the "wrong" sign, 
respectively. This should be helpful to differentiate from alternative charged pseu- 
doscalar decays leading to the same final states. Ultimately, the precise knowledge of the 
quantum effects on —>■ tb should provide the characteristic features necessary to pin 
down the supersymmetric nature of a heavy charged Higgs hypothetically discovered in 
the next generation of experiments at the Tevatron and at the LHC As we have shown, 
the mechanisms capable of producing a charged Higgs scalar in a hadron collider, e.g. tb 
and bt fusion, which can be greatly enhanced at large tan/5, are very sensitive to poten- 
tially large SUSY-QCD quantum effects. If these effects are eventually found, they could 
be the smoking gun needed to recognize that the produced in a hadron collider is. 



truly, a SUSY Higgs. After submiting this paper, we have noticed the work of Ref.[|18 
dealing with the same subject. 
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Figure Captions 

• Fig.l (a) SUSY-QCD Feynman diagrams, up to one-loop order, correcting the 
partial width of ^ tb; (b) Typical charged-Higgs production mechanism at 
hadron colliders. 
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Fig.2 (a) SUSY-QCD corrected T{H+ ^ tb) - eq.(^) - (in GeV) as a function 
of tan/?, compared to the corresponding tree-level width, Tq. Also shown are the 
partial widths of the alternative decays r+ z/^ and . The top- 

quark mass is mt = 175 GeV; (b) Comparative effects from the (finite) vertex and 
self-energies as a function of tan j3. The fixed parameters for (a) and (b) are given 
in the frame. 



Fig. 3 (a) Dependence of the relative SUSY-QCD correction Sg, eq.(0), upon the 
supersymmetric Higgs mixing mass parameter, /i; (b) Evolution of the SUSY-QCD 
correction in terms of the gluino mass for m^-^ = 65, 80, 100 GeV, tan P = 30. Rest 
of inputs as in Fig.2. 

Fig. 4 (a) Sg as a function of m^^; (b) Sg as a function of At,. In both cases tan/3 = 
10, 20, 30 and the remaining inputs are as in Fig.2 

Fig. 5 The 0{as) standard QCD correction, 6g, compared to the SUSY-QCD cor- 
rection, 6g, as a function of Mh+ and for two values of tan/3. Rest of inputs as in 
Fig.2. 
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